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Measurements of the suppression of the yield per nucleon and differential distributions of J/ij: 
production for 920 GeV/c protons incident on heavy nuclear targets have been made with 
broad coverage in pr and negative coverage in xf of produced meson. Production ratios of 
ip{2S) to J/ip and Xc to J/V^ have been measured with a high accuracy. The and 
production cross sections as well as to ratio have been obtained on one of the highest 
statistics available in proton nucleus experiments. 



1 Introduction 



The mechanisms by which charmonium and open charm states are produced in hadronic cohi- 
sions are not weh understood and are the focus of much interest. Although most of the results 
from previous experiments are qualitatively described by different theoretical models, there 
are some quantitative discrepancies that may be resolved with a better understanding of non- 
pert ubative effects in the production of heavy quarks. Taking into account that the production 
mechanisms can be identified by their strong kinematic dependences, it is crucial to have new 
measurements with high statistics and broad kinematic coverage to accomplish comprehensive 
descriptions of charm production in nuclei. In this context the HERA-B experiment provides 
new data on open and hidden charm production. 



2 The HERA-B Experiment 

HERA-B is a fixed target experiment at the HERA storage ring at DESY. Diverse hidden and 
open charm states are produced in proton-nucleus inelastic collisions. The pN energy in the 
center of mass corresponds to 41.6 GeV. 

The spectrometer emphasizes vertexing, tracking, particle identification and features a ded- 
icated multilevel J ftp trigger. 
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Figure 1: Layout of the HERA-B detector. 

The HERA-B target is a fixed target system consisting of two stations separated by 4 cm 
along the beam direction. Each station houses 4 target wires which can be moved independently 
into the beam halo. 

The side view of the detector is shown in Fig. The protons come from the left and 
interact in the target region. The target station is followed by the Vertex Detector System(VDS), 
which provides a precise reconstruction of primary and secondary vertices. The main tracking 
system(ITR+OTR) is placed behind the VDS. It uses different technologies, depending on the 
distance to the beam center. 

Particle identification is performed by a Ring Imaging Cherenkov Counter (RICH), an Elec- 
tromagnetic Calorimeter (EC AL) and a MUON detector (MUON). 

The J/^/J triggerll is initiated by di-lepton signatures from either MUON or ECAL systems. 
The information on regions of interest defined by pretriggers is used as seeds for the following 
search for tracks in the tracking system. Then the trigger performs the extrapolation of found 
tracks through the magnet and follows them through the VDS. In addition a vertex constraint 
is applied. Events with at least two fully reconstructed di-lepton candidates which come from a 
common vertex pass the trigger requirements. 

3 Data Set 2002/3 

The HERA-B collaboration collected data during 3 month data-taking period in 2002/2003. The 
data samples were obtained using carbon, titanium and tungsten wires. During that time 150 
million di-lepton triggered events were acquired. About 170,000 of J/'0 were reconstructed in 
muon channel and about 150,000 of J/^ were found in electron channel. The trigger performance 
has been significantly improved since the year 2000. For comparison, in 2000 the J/tp rate 
was 25-35 per hour and in 2002/2003 we already had 1000-1500 J/^ per hour. Significant 
improvements in the data acquisition system had been made as well, which resulted in achieving 
high recording rates, which in turn allowed accumulation of 210 million minimum bias events in 
just one week of running. The Minimum Bias trigger selected events satisfying a requirement on 
either a minimum number of hits in RICH or a minimum energy deposited in ECAL. It allowed 
to suppress events with no interaction without introducing any significant bias. 

4 Preliminary results. The Minimum Bias Data 

The following open charm decays are fully reconstructed in minimum bias data: 




Throughout this paper charge conjugate decays are also imphed. 

The invariant mass plots for and are shown in Fig. [21 and Fig. [HI respectively. The 
bump seen to the left of the in Fig. [HI is produced by reflection of not fully reconstructed 
charm decays, for instance K~7v^7v^. The shape of this reflection was obtained from MC 

generated cc events and was used in the fit as a component of the background in Ktv invariant 
mass plot. 




Figure 2: M{K7V7v). Figure 3: M{Kn). Figure 4: M{Knnsiow) - M{Kn). 

The mass difference between the mass of and the mass of is shown in Fig. 01 This 
quantity peaks close to the pion mass and has a very small width. The peak is narrow due to the 
fact that the errors in calculation of the invariant mass of three particles (in case of i?*^) and 
two particles (in case of D^) are highly correlated, and therefore cancel out to a large extent. 

The signal yields for the fully reconstructed ^ and decays are shown in Tabled 



Table 1: Signal yields for , and D* 



D+ K-Tr+7T+ 


^ K-7T+ 


D*+ 0*^%+ 


98 ±12 


189 ± 20 


43 ±8 



The sample of open charm collected by HERA-B constitutes one of the highest statistics 
available in proton beam experiments. Therefore the HERA-B will provide the precise measure- 
ment of D^^ and cross sections as well as ^ and ratios. 



Table 2: Results on and production cross sections. 





-0.1 <XF < 0.05 


full Xf 


(j{D^)lih/nucl 


11.5 ± ^-^stat i '^•'^sys 


30.2 ± A.^stat ± 5.8,^;, 


a{D^)lih/nucl 


2lA±?>.2stat^^-^sys 


56.3 ± 8.5stat ± Q-^sys 



The preliminary results on and cross sections are summarized in Tabled The xp 
range of this analysis extends from -0.1 to 0.05. Also shown are the total cross sections which 
were calculated by extrapolating over the full xp range. In the calculation of the ratio between 

and cross sections the systematic errors related to detector efficiencies are minimized and 
the dependence on the absolute luminosity determination is removed. This allows the accurate 
measurement of the cross section ratio. Result for the ^ ratio is 0.54 zb O.llgtat i 0.145^5. 

The comparison of our measurement for ^ ratio with those of other experiments is shown 
in Fig. ini For the ratio the experimental situation is unclear, and the uncertainties are large 



particularly for proton beam experiments. The line on the plot shows the expected theoretical 
value for the ratio of charged to neutral D mesons which is approximately equal to ^ and is the 
same for the proton and pion induced interactions ^. The big difference between the and 
production results from the fact that D mesons can be produced either directly or through 
feed-down from i?* decays. The prediction for the ratio is obtained by taking into account the 
well known branching ratios for the decays of i?* into D mesons. 



LEBC-EHS 
ACCMO^ I t 

- I M|lebc-ehs 

-ACCMb^ |-,gQ^\vA92 (|)E 653 



E769^ 



>E653 
HERA-B 



15 20 25 30 35 40 45 

EcMs(GeV) 







: DD{iih/nucl) 


HERA-B 
E653J * 

f 


LEBC-EHS 


E769 * 


LEBC-MPS 


; I 




' ACCMOR 

i 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

EBEAM(TeV) 



Figure 5: Ratio between and 
production cross sections. The line 
shows theoretical expectation. 



Figure 6: Comparison of our DD 
cross section with the results from pre- 
vious experiments. 



Our result on DD cross section is shown in Fig. El It is compatible with previous data 
and has smaller experimental uncertainty in comparison with the experiments at nearby proton 
beam energies. 



5 Preliminary results. The Di-lepton Data 

5.1 The Di-lepton invariant mass spectra 

The ji^ ii~ and e+e~ invariant mass spectra are shown in Fig. [71 and Fig. [HI respectively. Clear 
J and iIj{2S) resonance peaks are visible. Lower mass particle decays of p, cj, (j) into pair 
are also seen. 




Figure 7: invariant mass spectrum. Figure 8: e+e" invariant mass spectrum. 

Di-lepton trigger. Di-lepton trigger. 



5.2 J 1^ Differential Distributions 



The huge sample of J/'0 decays accumulated in 2002/2003 and the large acceptance covered 
by the HERA-B spectrometer provide a broad coverage in pt, up to 4.8 GeV. The differential 
distribution for J/'0 as a function of transverse momentum is shown in Fig. El Superimposed 
on the plot is a fit to the parameterization function: 



da 



A- 



1 + 



35 • TT • J)T 

256 ' <pt> 



Our results on the average transverse momentum as well as the results^^^ from E771 and 
E789 experiments are shown in Table 01 The results obtained in electron and muon channel are 
compatible with each other for each target material. The HERA-B measurement confirms the 
trend of increasing average pT with increasing atomic number and energy in the center of mass. 

The differential distribution of J/'0 versus Feynman xf \s plotted in Fig. E3 The points 
and errors are shown along with a fit to the parameterization function: 



da 

dxp 



A - {I - \xf\ 



The HERA-B differs from all other proton-nucleus collision experiments by its large accep- 
tance in the negative xp region. Our measurement for exponential parameter and the result^^^™^ 
from previous experiments are shown in Tabled Due to the ongoing acceptance estimations for 
different running periods, at the moment only a range for the slope is quoted. 
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Figure 9: J/'0 differential distribution versus pr- Figure 10: J/'0 differential distribution versus xf- 



Table 3: The HERA-B results on the average transverse momentum for J/V^. Comparison with data from E771 

and E789 experiments. 



Experiment 


Beam Energy 


Target 


Pt Range 


< Pt > {e^e ) 


<Pt> ) 


HERA-B 


920 GeV 


C 


PT < 4.8 GeV 


1.22 ± O.OUtot 


1.22 ± 0.01,tot 


HERA-B 


920 GeV 


W 


PT < 4.8 GeV 


1.29 ± O.OUtot 


1.30 ± 0.01,tot 


E771 


800 GeV 


Si 


PT < 3.5 GeV 




1.20 ±0.01 


E789 


800 GeV 


Au 


PT < 2.6 GeV 




1.29 ±0.009 



Table 4: The HERA-B results on the exponential parameter for J/'0 distribution in xf- Comparison with data 

from E771 and E789 experiments. 



Experiment 


Beam Energy 


Target 


xp Range 


c 


HERA-B 


920 GeV 


C,W 


-0.35 <XF < 0.15 


(5.0 - 6.5) ± O.Sstat 


E771 


800 GeV 


Si 


-0.05 <XF< 0.25 


6.54 ±0.23 


E789 


800 GeV 


Au 


-0.03 <XF < 0.13 


4.91 ±0.18 


E789 


800 GeV 


Cu 


0.30 <xf < 0.95 


5.21 ±0.04 



5.3 A-Dependence for J/t(j 

The suppression of charmonium production in heavy-ion colhsions is predicted to be a signature 
for the formation of the quark-gluon plasma. And it is therefore of particular interest to study 
these mechanisms in the absence of quark-gluon plasma. In this context the HERA-B presents 
the new measurement of A-Dependence for J ftp. 

Nuclear effects are usually considered in a power-law parameterization: 



The parameter a can be determined from measurements with two materials: 



1 



a = 



log{Aw/A^ 



c) 



log 



Nc Cw 



where N is the number of reconstructed J/^'s, L is the integrated luminosity, e - is the detector 
and trigger efficiency and the acceptance. To reduce the systematic effects only runs with C and 
W wires operating simultaneously were used. 

A-dependence exhibits strong kinematic dependences with Feynman xp and transverse mo- 
mentum pt of the produced J/'0. Therefore we plot a as a function of xp^ as shown in Fig. 
Since the ratio of luminosities is under investigation at the moment, the results were normalized 
to the E866 data^. 
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Figure 11: The A-dependence for J/^ip versus xf- Also shown are the results from E866 experiment. 



5.4 Production ratio of tl;(2S) to J/^Ij 

The production ratio of '0(25) to J/'0 is determined according to: 

^ a(V^(25)) ^ N{^jj{2S)) BrjJ/jj ^ l+r) ejj/j,) 
«V(2S) ^^j^^^ N{J/iP) ' Br{iP{2S) ^ 1+1-) ■ e{i>{2S)) 



The electron channel (see Fig. El gives 0.13 ± 0.02stat i^i the R^{2S) I'atio, which is in good 
agreement with the result from the muon channel. We also measure the ip{2S) cross section 
relative to the J/ip. The advantage of this method is that detector and trigger efficiencie s lar gely 
cancel. The value for the J/'0 cross section is based on E771 and E789 measurements'^^ and 
corresponds to a{J — 357 ± ^stat i 21 sysnh/nucl ^ . The resulting value for the ip{2S) cross 
section is a{t(j(2S)) = 46 =b 12 stat'i^b/nucl and is shown in Fig. along with the results from 
previous experiments. 



10^ 



10 

















2S) 


\ 






, , 1 , , , 1 , , , 1 , , , 1 , , , 1 , . , 1 , , , 1 , 


, , 1 , , , 1 , , , 



Z4 2,6 2.8 



3.2 3A 3.6 3.8 



70 
60 

50 
40 

50 
20 



E537 
_(W) 



i>E705 
i (Lf)... 



FNAL-B^k 
(C) 

FNAL-A 



FNAL-C 



(^^) !rE77l 



(Si) J. (C) 



e e InvariammassGeV/c 



, (GeV) 



Figure 12: The e e invariant mass distribu- 
tion. J/'0 and il^{2S) resonances are visible. 



Figure 13: The proton beam energy depen- 
dence of the ip{2S) production cross section. 



5.5 Fraction of J/tl^ produced via Xc decays 

The production ratio of Xc to J ftp can help to discriminate between different theoretical models. 
The Xc is reconstructed in the decay Xc J /i^l ~^ Due to the small branching ratio of 

XcO J/'^^li only Xci and Xc2 contribute to the observed Xc signal. The energy resolution in 
the HERA-B is insufficient to distinguish Xci and Xc2 states. The fraction of J/^ produced via 
radiative Xc decays is determined from: 

The published result based on data from an earlier HERA-B run with limited statistics for 
Xc is Rxc — 0-32 ± O.OQstat i 0.045^5^^, which is compatible with most of the previous data and 
the NRQCD prediction. 

In the 2002/2003 HERA run, HERA-B collected a Xc sample which significantly exceeds the 
statistics obtained by previous experiments. The mass difference distribution AM = M(/+/~7) — 
M{l^l~) is shown in Fig. El An excess of events with respect to combinatorial background 
determines the number of Xc candidates. The same plot, but after background subtraction is 
shown in Fig.QSl The shape of the background is taken by combining J/'0 and photon candidates 
from different events. The plot corresponds to 15% of statistics available in muon channel. We 
expect about 15,000 Xc foi" the full data set combined from electron and muon samples. The 
preliminary result for R^^ from the 2002/2003 data is 0.21 ± 0.05stat- 

6 Summary 

HERA-B has finished data-taking and now focuses on analysis of large accumulated data sam- 
ples: 150 million of events with di-lepton trigger and 210 million of events with minimum bias 




Figure 14: The AM = M(l^rj) - M{l^r) Figure 15: The Xc signal after background 

distribution. subtraction. 



trigger. These data were obtained on a variety of nuclear targets. About 300,000 J/^'s decay- 
ing into electron and muon pairs were reconstructed in the di-lepton sample. HERA-B features 
broad coverage in pt and explores the negative xp region for the first time. Differential distribu- 
tions for J/ip as functions of pt and xp are shown. New data on J/ip suppression are presented. 
Results on the production ratios of i^(2S) to J/^ and Xc to J/tl) are also shown. Clear signals 
of and mesons were obtained in the minimum bias data. Results on and 

production cross sections and the to ratio are discussed. 

Lots of physics analyses are ongoing and we are looking forward to new results. 
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